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Abstract 

We investigate the thermodynamics of a four-dimensional charged black hole in 
a finite cavity in asymptotically flat and asymptotically de Sitter space. In each 
case, we find a Hawking-Page-like phase transition between a black hole and a 
thermal gas very much like the known transition in asymptotically anti-de Sitter 
space. For a "supercooled" black hole — a thermodynamically unstable black 
hole below the critical temperature for the Hawking-Page phase transition — 
the phase diagram has a line of first-order phase transitions that terminates in 
a second order point. For the asymptotically flat case, we calculate the critical 
exponents at the second order phase transition and find that they exactly match 
the known results for a charged black hole in anti-de Sitter space. We find strong 
evidence for similar phase transitions for the de Sitter black hole as well. Thus 
many of the thermodynamic features of charged anti-de Sitter black holes do 
not really depend on asymptotically anti-de Sitter boundary conditions; the 
thermodynamics of charged black holes is surprisingly universal. 
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1 Introduction 



It has been known for some time that a black hole in asymptotically anti-de Sitter 
space can undergo a phase transition, "evaporating" at a critical temperature into a 
hot gas pQ. This Hawking-Page phase transition has recently garnered new attention 
with the realization that it can be used to test the AdS/CFT correspondence by 
identifying a similar structure in a dual conformal field theory [2 . A dual description 
may exist not only for this particular transition, but for a number of finer details of 
the phase structure of the asymptotically AdS black hole as well [HlUj . The purpose of 
this paper is to show that this phase structure is not unique to the asymptotically AdS 
black hole, but occurs universally for suitably stabilized black holes in asymptotically 
flat and asymptotically de Sitter space. 

A black hole in anti-de Sitter space can be thermodynamically stable. In asymp- 
totically flat or de Sitter space, this is no longer the case: an isolated black hole 
radiates away its energy in the form of Hawking radiation. To understand equi- 
librium black hole thermodynamics, one must therefore work with ensembles that 
include not just the black hole, but also its environment 0. As self-gravitating sys- 
tems are spatially inhomogeneous, any specification of such ensembles requires not 
just thermodynamic quantities of interest, but also the place at which they take the 
specified values. Having set up an appropriate thermodynamic ensemble, one can 
then proceed to ask interesting questions regarding stability and phase structure. 

Such an analysis was done for the case of charged black holes in asymptotically 
flat space by Braden, Brown, Whiting, and York [jy. Following a brief summary of 
their analysis, we will build on those results to study the phase structure of charged 
black holes. We will also apply a similar method to charged black holes in de Sitter 
space. This work grew out of an effort to understand the effect of varying "con- 
stants" on black hole thermodynamics 0, but it may have broader significance. In 
particular, the universality of the phase structure may have interesting implications 
for holography in asymptotically flat and de Sitter space. 

The Euclidean action [H] for metric g and electromagnetic field over a region 
M with boundary dM takes the form 

J = --L/" d A x^-g[(R-2k)-F 2 ] + ^- f d 3 x^(K-K ), (1.1) 

1D7T Jm OTT JdM 

where the cosmological constant A is zero for flat space and positive for de Sitter 
space. The geometric and electromagnetic data specified at the boundary dM fixes 
for us the thermodynamic data for this gravitational system thought of as a statistical 
ensemble. We take a spherically symmetric ansatz for the metric of a black hole 
spacetime, 

ds 2 = f 2 dr 2 + a 2 dy 2 + rW, (1.2) 

where /, a and r are functions of the radial coordinate y £ [0, 1]. The reduced action 
may be calculated from 1)1. 1J) provided certain conditions are taken into account: 
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1. The horizon of the black hole is located at y = (i.e. r(0) = r+), and we 
restrict to the class of metrics that are "regular." In other words, the geometry 
of the spacetime near the origin of the y-r plane looks like the plane M 2 , so the 
horizon is nondegenerate. 

2. The boundary at y = 1 has a topology S 1 x S 2 with the two-sphere having 
an area 4-7rr^. Heat can flow in either direction through this boundary in such 
a manner as to keep the temperature (3~ l fixed. The inverse temperature is 
simply the proper length of the circle S 1 of the boundary: j3 = 27r/(l) (where 
we take r to have period 2tt). 

We also need to specify appropriate electromagnetic data at the boundary y = 1. For 
the grand canonical ensemble, we specify the gauge potential A T {\) (more precisely 
the difference of potential between y = and y = 1). For the canonical ensemble, we 
must add a boundary term to the usual Maxwell action and hold the electromagnetic 
field fixed at the boundary. This is equivalent to keeping the total charge inside the 
cavity fixed. 

Using this analysis, we first demonstrate in section 2 that black holes in asymp- 
totically flat space show evidence of a Hawking-Page phase transition to a thermal 
gas. While the transition temperature depends on the charge, the variation is sur- 
prisingly small. In the absence of a well-understood comparison action for a charged 
gas, conclusions from (fixed-charge) canonical ensemble are somewhat ambiguous; we 
therefore also check for this transition in the grand canonical ensemble, and confirm 
its existence. 

We next investigate the unstable "supercooled" region of the phase diagram, and 
show that black holes in this region, like those in anti-de Sitter space [130], exhibit 
a first order phase transition: for a given (small) value of charge, the entropy of the 
system changes discontinuously as the temperature is varied. As the charge increases, 
the "height" of this discontinuity decreases, and the line of first order phase transitions 
terminates at a critical point. We show that this critical point is the location of a 
second order phase transition, and we calculate the critical exponents associated with 
it. Remarkably, the critical exponents are the same as those found by Chamblin et 
al. |3]li| for charged black holes in anti-de Sitter space, hinting that these two systems 
belong to the same universality class. 

In section 3, we apply the reduced action technique to charged black holes in 
de Sitter space. We again find evidence for a Hawking-Page phase transition. For 
sufficiently small cosmological constant, the supercooled region again displays a line 
of first order phase transitions terminating at a critical point. This picture no longer 
remains valid as A increases. In particular, it seems that for large values of A, there is 
no second order transition: the location of the would-be transition is pushed beyond 
the allowed boundaries of the cavity. 

The reduced action for the de Sitter case is more complicated, and we cannot find 
the critical exponents of the second order phase transition (when it exists) analyti- 
cally. We can study two special cases: the uncharged black hole for any value of A and 
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the charged black hole for small A. We again find the same structure, with critical 
exponents that are unchanged by the presence of a small cosmological constant. 
Section 4 is devoted to discussion our results. 




2 Charged black hole in asymptotically flat space 

Following Brown et al. [§], we first insert the ansatz (|1.2|) into the action 1)1.1(1 . and 
find that the Hamiltonian constraint can be solved, yielding 

V(r)=f-V = l-- + 4- (2-1) 

For a nondegenerate horizon, C = r+ + e 2 /r+, with e 2 < r\. The reduced action for 
the canonical ensemble then takes the form 

Ic = Pr B [l-,l(l-^)(l--£-])-,rl. (2.2) 

Extremizing the action with respect to r + gives us 

^^-^-^(l-^) 1 ' 2 , (2.3, 

which should be viewed as an equation for the unknown r + in terms of the fixed 
charge, boundary radius, and temperature of the canonical ensemble. Each extremum 
gives a saddle point contribution to the Euclidean path integral for the canonical 
partition function Z[/3], with a Helmholtz free energy and entropy 

F = ~\nZ[0\ = ^I c %r + {(5)\ 

S = -(/3^-l)lnZ[/?]=^r 2 . (2.4) 

Equation (|2.3|) can be rewritten as the seventh order algebraic equation 

x 5 (x-l)(x-q 2 )+b 2 (x 2 -q 2 ) 2 = (2.5) 

with x = r + /rg,q = e/rg and b = ft/Anr^. We are interested only on those solutions 
that are real positive and satisfy q 2 < x 2 . Of these positive solutions, we initially 
choose the one that minimizes the free energy F[/3], since this is the thermodynami- 
cally stable state in the canonical ensemble. 
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Figure 1: Free energy of a charged black hole in asymptotically flat space as a function 
of inverse temperature b and charge q. 

2.1 The Hawking-Page transition 

The plot of free energy as a function of charge q and inverse temperature (3 is shown 
in Figure ^ For reference, we have also plotted the plane F = 0. In the original 
paper of Hawking and Page pQ, the phase transition to a hot gas of particles was 
determined by comparing the black hole free energy with that of a real thermal gas, 
but subsequent work (for example, has often used "hot empty AdS" as a reference. 
This alternative is not available for a charged black hole in the canonical ensemble, 
except at q = 0, since there is no "hot empty charged space" ; one cannot avoid the 
intricacies of a real charged gas. 

Although there has been a bit of work on the thermodynamics of self-gravitating 
gases of particles [^JEH^D; ^ e surj j ec1, remains rather poorly understood. For low 
temperature, Refs. PllOpilj confirm the natural guess that the free energy is approx- 
imately Nm, where N is the number of gas particles of mass m and charge qo, while 
the charge is e = Nq^.* For real elementary particles, q$^> mm. geometrized units — 
for an electron, for instance, qo/m ~ 10 21 — so FCe, even for a near-extremal black 
hole. An approximation F = for a thermal gas thus seems reasonable, and the 
intersection of the F = plane with the free energy in Figure ^ should give the rough 
location of a Hawking-Page-type phase transition. It is straightforward to check that 
above the transition temperature, the black hole is thermodynamically stable, i.e., 
Cy = —b(dS/db) is positive everywhere. Figure |2J for example, plots the entropy 

*For Newtonian results, see also 
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Figure 2: Entropy of a charged black hole in asymptotically flat space as a function 
of inverse temperature b and charge q. 

against q and (3; the sign of Cy is determined by the slope. 

In the absence of a better understanding of the thermodynamics of hot self- 
gravitating gases, this argument is admittedly not completely convincing. The "com- 
peting" configurations are, after all, not just dilute gases, but stars. To check more 
carefully for a Hawking-Page transition, we have therefore followed the suggestion of 
Ref. JH] and looked at the grand canonical ensemble, in which the boundary poten- 
tial (j) is fixed instead of the charge. This ensemble has the nice feature that it does 
allow a comparison with "hot empty space": while an empty cavity cannot contain 
charge, its walls can be held at a constant potential. 

Figure |H1 shows the Gibbs free energy of the grand canonical ensemble, with the 
plane G = superimposed. The intersection shows the location of the Hawking- 
Page phase transition to hot empty spaced The transition depends only weakly on 
the potential. There is one important difference between the charged and uncharged 
cases, however. For an uncharged black hole, there is a further critical temperature 
beyond which the action has no extremum, even with positive free energy. Thus 
even a "supercooled" black hole cannot exist below this temperature. For a charged 
black hole, no such temperature exists; a locally stable black hole can exist for any 
value of /3, essentially because a near-extremal black hole can have an arbitrarily low 
temperature. 

^The Gibbs free energy becomes positive again at very high temperatures. This occurs when the 
thermodynamically stable black hole becomes larger than the cavity radius tb, leaving only an 
unstable (negative specific heat) configuration. 
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Figure 3: Gibbs free energy of a charged black hole in asymptotically flat space as a 
function of inverse temperature b and potential (p. 



2.2 Supercooled black holes and critical exponents 

Let us now return to the canonical ensemble, and consider the portion of Figure ^ 
for which the black hole free energy is positive. While a black hole in this region is 
globally unstable, it is locally stable — the specific heat is positive, as is evident from 
Figure 12 One can therefore produce a "supercooled" black hole in this region, and, 
indeed, the analysis of fluctuations for uncharged black holes suggests that such a 
state may have a very long lifetime |141ll5j . Figure ^ now suggests that there is a 
line of first order phase transitions in this region that terminates at a critical point 
in the (q, /3)-plane. We will show that this critical point is the location of a second 
order phase transition. 

The simplest way to see this is to rewrite ()2.3|) as 

l(M) = !»|W!. (2 . 6 ) 

1 — q A /x A 

In general, b(x, q) has may have extrema (both maxima and minima) as a function 
of x, with locations 

db _ x(<? 4 (5 - 6x) + (2 - 3x)x 3 + 2q 2 x{-3 + 3x + x 2 )) _ 

fa ~ 2{q 2 - x 2 ) 2 v / T^V 1 - I 2 /* ~ 

The numerator is a polynomial in x (with functions of q as coefficients), with a 



discriminant 



A( 9 ) = ^ 4 (l-^) 4 (l-18g 2 + g 4 ). (2.8 
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This discriminant determines the condition for the maximum and the minimum to 
coincide: they do so for A(q) = 0, i.e., for q c = —2 + y/5. The corresponding critical 
values for the horizon radius and temperature are x c = 5 - 2^5 and b c = §(17 + 
ZS/Vly 1 / 2 . Thus (x c ,b c ,q c ) ~ (0.528,0.429,0.236). Expanding in the neighborhood 
of the critical point, we obtain 

b-b c = --(425 + 195\/5) 1/2 (x - x c f + ■ ■ • ~ -3.644(x - x c f + ■■■ . (2.9) 
8 

To see the dependence on x near q c , we rewrite (|2.3j) as 

^ _ (26V - x 5 + x 6 - (1 - x)xV46 2 + x 4 ) 1 / 2 

6) - (2.10) 

and expand in the neighborhood of (x c , 6 C ): 

q-q c = ^j(105 + 47\/5)(x - x c ) 3 + • • ■ ~ 5.252(x - x c ) 3 + • • • . (2.11) 

Inverting these relations and writing in terms of the entropy S, we find 

S-S c = 1.226(r - rj 1 / 3 + • • • , (2.12) 
S-S c = 41.908(g - q c ) 1/3 + ■ ■ ■ , (2.13) 

where r = rskT,T c = l/b c and 5c = 7rx 2 . The specific heat = T(dS/dT) = 
T(dS/dr) is 

C v ~ 0.409(r - r c )" 2/3 (2.14) 

which means that the critical exponent a is —2/3. Remarkably, this is also the critical 
behavior of the charged black hole in asymptotically anti-de Sitter space |3*lHlll8|. 
There, as in the asymptotically flat case discussed here, a line of first order phase 
transitions in the "supercooled" region ends in a second order transition, and although 
the AdS partition function has a completely different functional form, the critical 
exponents are identical. 

The case of uncharged black hole is even simpler. Equation (|2.5j) reduces to the 
cubic x 3 —x 2 +b 2 = 0, with critical values (x c ,b c ) = (2/3,2/(3\/3)). This critical point 
does not signal a second order phase transition in contrast to the case of even 
arbitrarily small charge, the action for the uncharged black hole has no extremum 
beyond this point. As discussed at the end of section |2~T1 this point instead describes 
a transition between even a supercooled black hole and a thermal gas. Near this 
point, 

S-S c ~ (r-r,) 1 / 2 (2.15) 

indicating that the specific heat diverges like C v ~ (r — r c ) -1 / 2 . This is that same 
critical behavior seen in non-spinning L>3-branes |17|I19|. 
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3 Charged black holes in de Sitter space 



Conceptually, the method of analyzing the thermodynamic stability of black holes 
in de Sitter space is the same as in flat space. Lorentzian de Sitter space has in- 
finite spatial extent, but in static coordinates there is a cosmological horizon. The 
corresponding Euclidean section of the metric has finite volume, with topology of a 
four-sphere with a radius equal to the cosmological horizon. 

For a charged black hole in de Sitter space, the spacetime has, generically, three 
horizons: an inner and an outer horizon for the black hole, and the cosmological 
horizon. Any discussion of equilibrium black hole thermodynamics requires one to 
specify thermodynamic quantities at a boundary that is inside the cosmological hori- 
zon. In this case, the entropy of the equilibrium configuration depends on three 
parameters: temperature, charge, and the cosmological constant. Again using the 
spherically ansatz (|1.2|) in the action and solving the Hamiltonian constraint, we find 

v (r) =(!iy=i + 4-£-£. ( 3.D 

\a J r L 6 r 

The integration constant C is often identified as the mass parameter 2M. 

Our first task is to determine the allowed values of the parameters. For the 
path integral to be well-defined, we need to restrict to those values of e, M, A and 
r for which V(r) is nonnegative. Event horizons are located at those values of r 
for which V(r) = 0. To locate these, it is useful to define dimensionless quantities 
p = ry/A, g = evA and p = My/A, so 

V = --^ Q(p; H, g 2 ) = (p 4 - 3p 2 + 6pp - 3g 2 ) (3.2) 

The quartic polynomial Q has four real roots: one negative root, and three pos- 
itive roots that we denote, in ascending order, p-(p,g 2 ), p + (p,g 2 ), and p c (p,g 2 ). 
The region outside the outer horizon of the black hole is thus p+(p, g 2 )/y/A < r < 
p c (p, g 2 )/y/A. The condition for two (or more) of these roots to coincide JB] is given 
by the vanishing of the discriminant of Q(p, \x, g 2 ): 

A(p, g 2 ) = (1 - Ag 2 f - (1 + Ug 2 - ISp 2 ) 2 = 0. (3.3) 

The the allowed values of mass and charge parameter in the (p, g)-plane are thus 
given by the interior of the region bounded by the curves 



1 + 12^-^/(1 -4g2)3 

Vmin{g) = V 18 ' ( 3 - 4 ) 



Hmax(g) = J 1 + U9 ^ (3.5) 

which intersect at the point (p,g) = (y/2/3, 1/2). The curve p m in{g) corresponds 
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Figure 4: Allowed values of mass fi and charge g. 



to r_ = r_|_, i.e., coincident inner and outer horizons, while the curve Pmax{d) cor- 
responds to r + = r c . These two conditions can be combined into a simple form 
g 2 = P 2 (l- P 2 ). 

The de Sitter black hole action is thus a function of four variables, p, g, b = /3\/A, 
and the box size R = r B \^A. To locate a saddle point, we must find the roots of 
V(r) = 0, choose the middle positive root, substitute it into the action, and minimize 
it, subject to the condition that p and g lie in the interior of the region depicted in 
Fig.H 

To analyze the stability of the black hole, it is convenient to revert to the variables 
r_|_ and e rather than their "normalized" versions p and g. We will also stay away the 
extremal black holes corresponding to the curves p m in{g) and p ma x(g), by requiring 
that g 2 < p\{1 — p\) and g 2 < p 2 (l — p 2 ). Requiring that the r-r plane look like M 2 
near r = r + fixes C to be r + + e 2 /r + — Ar^/3, which puts V(r) in the form 



V(r) 



l- r -± 



1 



l /r H 



A. , o n 

— (r + rr + + 
3 



(3.6) 



Thermodynamic quantities are specified at a boundary r = r#, and we require that 
r _B < r c- For any given value of r + and e, the maximum allowed value of is given 
by the larger root of the equation g 2 = p 2 (l — p 2 ). To summarize, for any given value 
of x = r + /rs, the allowed values of charge q = e/rs are given by q 2 < x 2 (l — R 2 x 2 ) 
and q 2 < 1 - R 2 . 

Subject to the above conditions, the reduced action can be shown to be 

2 



Ic 



pr B [l-^/V(r B )]-Trr' 




A 



r B r + 



{r B + r B r + + r\ 



(3.7) 
vrr 2 .(3.8) 
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(a) R = 0.5 (b) R = 0.9 

Figure 5: Helmholtz free energy of the black hole for two values of R. 

Extremizing that action with respect to x and solving for b, we obtain 



4vrr + ^(l-^)(l-^-f(r|+r + r B+ ^ 

^ ' 1 " 4~ B ~ f (4 + r + r B + r* ) + (l - g) (f (r| + 2r + r B ) - |) ' ^ 
This is a ninth order algebraic equation for r + , which may be written as 

X 5 (l - X) ( X - q 2 - R2x ( 1 + X + x2 ) \ _ b 2 {g 2 _ x 2 + R 2 X ^2 = _ (3 1Q) 



Given values for b, q and R, (|3.1U|) can be solved numerically for x. In Figure |S1 
we plot the free energy as a function of charge q and inverse temperature b for two 
values of the cavity radius R. As before, we also show the F = plane, which should 
give an approximate locus for the Hawking-Page phase transition. 

As in the asymptotically flat case, we have checked this behavior in the grand 
canonical ensemble. Figure El shows the Gibbs free energy and the G = plane, 
and confirms the existence of a phase transition. Note that a Hawking-Page phase 
transition occurs for all values of A. While the presence of a cosmological constant 
affects the transition temperature, it does not change the qualitative behavior. 

For small values of R, it is clear that the charged black hole in de Sitter space 
behaves very much like the corresponding black hole in flat space. There is again a 
line of first order phase transitions in the supercooled region that ends at a critical 
point. For large values of R, the critical point disappears: the place at which it would 
occur is essentially pushed outside the cosmological horizon. 

In general, we have been unable to find the nature of the phase transition at 
the critical point. We can, however, do so for for two important special cases: the 
uncharged black hole (e = 0) and the case of small cosmological constant. For e = 0, 
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(a) R = 0.3 



(b) R = 0.5 



Figure 6: Gibbs free energy of the black hole for two values of R. 



the function b takes the form 

b(x, R) = *M-M-J^+*±*!m. (3.1!) 

The positive root of the equation 1 — A(r 2 + r + r + r+)/3 = gives us the location of 
the cosmological horizon r c . Requiring that r + < r c tells us that r\ < 1/A. Similarly, 
r + < tb < r c implies that R = r£\/A < 1. 

The free energy can be plotted as a function of b and R (Figure and clearly 
shows that black holes do not exist above a critical value b c of b, for any fixed value 
of R. This critical value is analogous to the value discussed at the end of section 
12.11 it is a limiting temperature beyond which not even a supercooled black hole can 
exist. For any given value of R, the entropy of the system behaves as 5 — S C (R) ~ 
(b — b c (R)) a . By choosing various values of R, one can check that a = 1/2 and that 
Cv ~ ( T — T c]~ 1 ^' 2 ■ The critical behavior thus is very much like that of the uncharged 
black hole in flat space discussed in section l2~2l 

For the case of small A, we can expand b in powers of R 2 and look for critical 
points in the supercooled regime. Both b c and q c receive correction to order R 2 , while 
x c remains the same to this order. The critical exponents also remain the same. 



4 Discussion 

The thermodynamics of black holes is interesting in its own right, and also because 
of its connections to holography and the AdS/CFT correspondence. Largely inspired 
by the latter connection, past investigations have found a rich phase structure for 
asymptotically anti-de Sitter black holes. 

Our main result is that this phase structure does not, in fact, require asymp- 
totically anti-de Sitter boundary conditions, but appears as well for black holes in 
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0.5 

Figure 7: Free energy for the uncharged de Sitter black hole as a function of inverse 
temperature b and cosmological constant R = y/Xrs- 

cavities in asymptotically flat and asymptotically de Sitter spaces. In particular, both 
Hawking-Page-like phase transitions and transitions between small and large black 
holes occur as we vary the charge and temperature. For "supercooled" black holes 
in flat space, and black holes in de Sitter space with small A, we found a line of 
first order phase transitions ending at a second order transition, and we were able to 
identify the critical exponents at the second order point. Remarkably, these were the 
same as those for asymptotically anti-de Sitter black holes, despite a very different 
algebraic structure of the partition function. 

For anti-de Sitter space, the phase structure of charged black holes has served as 
a valuable test of the AdS/CFT correspondence. We leave open the very interesting 
question of whether the corresponding results for asymptotically flat and asymptoti- 
cally de Sitter black holes might cast light on proposals for finite- volume holography 
in this more general setting. 
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